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Abstract 

Functions satisfying a defective renewal equation arise commonly in 
applied probability models. Usually these functions don't admit a explicit 
expression. In this work we consider to approximate them by means of 
a gamma-type operator given in terms of the Laplace transform of the 
initial function. We investigate which conditions on the initial parameters 
of the renewal equation give optimal order of uniform convergence in the 
approximation. We apply our results to ruin probability in the classical 
risk model, paying special attention to mixtures of gamma claim amounts. 
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1 Introduction 



For a given interval / £ R, let C{I) be the class of continuous functions g : 
/ — > R. The aim of this paper is to study the approximation of a function 
g G C([0,cxd)) in terms of its Laplace transform. To this end, we assume that 
\g{u)\ — 0(e^") as li ^> oo, for some 7 > 0. Then, the Laplace transform of g 

/>oo 

g{t) / e-*"9(u)d«, t>j (1) 
Jo 

is well defined and infinitely differentiable. Moreover, we can approximate g in 
terms of the derivatives of its Laplace transform. From now on, for a given func- 
tion g, g^"-^ will denote its n-th derivative (3^°^ :— g). We define the following 
operator 

L;g{u) = , u > 0, t > 7 (2) 

i {[tu\ + 1) 

where [u] indicates the largest integer less than or equal to u and r(-) is the 
gamma function. The approximation properties of ([2]) , as t 00, can be studied 
taking into account that the previous formula admits the following representa- 
tion (cf. 2, example (c), p. 92]). Let {S{t), t > 0) be a collection of random 
variables such that S{0) = and, for each t > 0, S{t) is a gamma r(i, 1) random 
variable. Recall that a gamma r(a,/3) random variable has density given by 

/a,M") - f^/3"u"-ie-^", M > /? > 0, a > 0. (3) 

By differentiation under the integral sign in ((T)) it can be seen that 

L;g{u) = Eg(^M±^], u>0,t>j (4) 



t 

A modification of the operator defined in ^ was used in [9] to approximate the 
distribution function Fx of a nonnegative random variable X by means of its 
Laplace-Stieltjes transform, using that when g — Fx can be rewritten in 
terms of a well-known inversion formula for this transform (in the beginning of 
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Section 3.2. we explain this connection). By considering the inversion formula 
^ using the Laplace transform instead of the Laplace-Stieltjes transform, we 
can therefore widen the class of functions under consideration, using, at the 
same time the general convergence results given in [5] . The connection of both 
inversion formulas through the same operator was considered also in [2]. Also, 
it is interesting to point out that similar approximation formulas involving the 
Laplace transform, have been used in the literature in order to obtain results 
concerning characterizations of life distributions in reliability and shape prop- 
erties of renewal functions (see |5J and the references therein). Moreover, in a 
recent paper (cf. i6j, we can find interesting numerical comparisons for dif- 
ferent inversion formulas involving the Laplace-Stieltjes transform of measures 
concentrated on the positive semiaxis. It should be mentioned that applied 
to distribution functions, is the so-called Widder's formula in [B] (in Section 3.2. 
we give more details). 

Remark 1.1 The well-known Post Widder's inversion formula for Laplace trans- 
forms, that is, 



admits a similar probabilistic interpretation (cf. Feller p.233]j, as we have 



Observe the main differences of the inversion formula ([2]) with Post- Widder 
inversion. In the first one, for t fixed, the order of differentiation increases with 
u, whereas the point at which the Laplace transform is applied remains fixed. In 
Post- Widder inversion, for t fixed, the order of differentiation is fixed, whereas 
the point at which the Laplace transform is applied varies with u. 

Our aim in this paper is to construct, as in [5] an accelerated approximation 
to ^ and apply it to functions satisfying a defective renewal equation of the 




(5) 




(6) 
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form: 

nu 

m{u) — (f> / m{u — y)dF{y)dy + v{u), u>0, (7) 

"'0 

in which F is the distribution function of a nonnegative random variable with 
F(0) = 0, is a parameter such that < < 1 and v : [0, oo) ~> R is a 
locahy bounded function. Many quantities of interest in applied probability 
(ruin probability in insurance theory for instance) satisfy a defective renewal 
equation. For specific references, along with the properties we are going to 
mention below, see for instance [Hi p. 152]. It is known that there is a unique 
locally bounded solution of ([7]). In fact, let F*" the n-th convolution of F with 
itself (F° being the point mass at 0). Define G := E^r=o(l " ^)'/'"^*"- Then, 
the above mentioned solution to (O is given by 

"^(") = T^/ viu~y)dG{y)+v{u), u>0. (8) 

^ ^ J{0,u] 

Only in very specific situations one can find an explicit solution to m using (|5]). 
From now on, we will assume that F is absolutely continuous, having density 
/, and therefore, ^ becomes 

m{u) = (j) / m{u — y)f{y)dy + v{u), u>0, (9) 
Jo 

Note that in this case the Laplace transform of m can be written, using ([9]) as 

rh{t) = ""^^l (10) 
1-0/(0 

and we can build the approximation i^m defined in ([2]), in cases that m cannot 
be computed in an explicit way. 

The paper is organized as follows. In the next Section we introduce a modi- 
fication of the operator defined in ^ improving the rate of approximation and 
apply it to renewal functions, investigating conditions under which the rate of 
uniform convergence is optimal. In Section 3 we consider a particular applica- 
tion of or our results in the context of ruin theory in insurance risk models. 



3 



2 The accelerated approximation. Application 
to renewal functions 



In order to improve the rate of approximation given by ([2]), we will consider, 
in a similar way as in ^ the following accelerated approximation for a given 
function g G C([0, oo)) 



I .9(0), iffc = 0; 

and for u > such that u ^ k/t, fc = 1, 2, . . . , 

Mf = {tu - N)Mf + (N + 1 - ^^^)Mf / . (12) 



improves the initial order of convergence of L^g (at most 1/t) to for 
suitable functions. In particular, in 9 a class of functions was studied under 
which this order of convergence holds uniformly. We now introduce this class 
of functions. 

From now on we will denote by C"([0, oo)) the subclass of functions in 
C([0, oo)) having continuous n-th derivative on (0, oo). Also, for a given function 
g : / — > R, \\g\\ will denote its uniform norm, that is 

hW :=sup|5(m)|. 

For a given subset A C /, we will use the notation \\g\\A '■= sup„g^ l5(w)|. We 
introduce the following class of functions: 

Pi := {g G C4([0, oo)) : \\g"{u)\\ < oo and \\u^g'^{u)\\ < oo}. (13) 



Remark 2.2 For g e Vi we have that ||ug"'(u)|| < Wu^ g'''' {u)\\ < oo (cf. [H 
p.571]j. The quantity \\ug"'{u)\\ will also appear in our error bounds. 
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The following result will play an important role from now on. 

Theorem 2.1 ([9], p. 571) Let g G Vi, with Vi as defined in and let 
Mj'^'g, t> be as defined in (ITT])-([n]). We have 



Remark 2.3 It would be interesting to compare the accelerated inversion for- 
mula Ml g with a similar procedure for the Post-Widder inversion formula 
given by 

Gf^giu) ■.^2W2tg{u)-Wtg{u). 

Previous expression is the classical Stehfest enhancement of order two for the 
PostWidder formula. A numerical comparison example between both accelerated 
inversion formulas was given in ^9, Example 2.1. p. 564]. The test function 
considered was g{u) = 1 — (1 — p)e^P",M > 0, and the parameters taken were 
p = Q.\ and t = 5. The advantage of this function is that both approximations 
can be easily computed. The numerical performance of each method for different 

[21 

values of u can be seen in 9, Table 2, p. 567], being Mf' g more accurate, 
specially for big values of u. Roughly speaking, the better accuracy of g{u) 
can be explained by the fact that the variability of the underlying random variable 
defining L^g(u) (recall (|4]) ) has, for fixed t and u — oo, less variability than the 
one defining the Post-Widder operator (recall On the other hand, the use 

of instead of L^g, improves substantially the numerical performance of 

the approximation for values u close to the origin (as noticed in [6], Widder's 
formula has less precision for values u close to the origin). See [9, Example 2.1. 
p. 564] for a more detailed discussion. 

Our aim is to consider the renewal function given in ^ in order to obtain 
conditions on F and v such that m e Pi , with 2?i as defined In this case, 
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by Theorem 12.11 Af| 'm, as defined in ((Tl]) - p2)) has order of convergence 1/i^. 
To this end, we need suitable expressions for the derivatives of m. From now 
on, we will denote by Cq ([0, oo)) to the subclass of functions in C"([0, oo)) such 
that lim q^'^Ht) exists for all fc = 0, 1, . . . n and is finite. 

First of all we state a technical lemma in order to justify differentiation under 
the integral sign in expressions similar to ([9]) . This result will be systematically 
used along the paper. 

Lemma 2.1 Let c : (0, oo) — R &e a function satisfying 

c(u) — / a{u — z)h{z)dz, it > 0, 
Jo 

in which a £ Cq([0,oo)) and b G C((0, oo)) is such that \b{z)\dz < oo for all 
u > 0. Then c is difjerentiable for all u > and verifies 

c'{u)= I a' (u - z)b{z)dz + a{0)b{u) 
Jo 



Proof. Consider u > 0, let < S < u and let < \h\ < S. To show the assertion 

we will check that 

1 / /-ti+'i \ 
lim— / a{u + h — z)b(z)dz — / a(u — z)b(z)dz \ = / a' (u—z)b(z)dz+a{0)b(u) 

h^o h yjo Jo J Jo 

Define a{u) — a'{0)u + a(0) if —S < u < (in this way a{u) is continuously 
differentiable on [—5, oo)). Observe firstly that we can write 



1 
h 



a{u + h — z)b{z)dz — / a{u — z)b{z)dz 
Jo ) 

— {a{u + h — z) — a{u — z))b{z)dz ^ / a{u + h — z)b{z)dz (14) 

on- nJu 

To deal with the first term, note that when h < and z varies between and 

u, then u + h — z varies between h and u + h, thus taking negative values. 

That is why we defined a on [—3, oo). Notwithstanding, as a is differentiable on 



this interval, and ~5 < h, we have always \h~^{a{u + /i — z) — a{u — z))b{z)\ < 
W'^' \\[-s,u+s]\b{z)\. As the right-hand side is an integrable function, we can apply 
dominated convergence theorem to write 

lim / — {a{u + h — z) — a{u — z)) b{z)dz ~ / a {u — z)b{z)dz (15) 
h^Q Jq h Jo 

As for the second term in (jl4[) . we can write 

— / a{u+h— z)b{z)dz — — {a{u+h— z)—a{u— z))b{z)dz+— / a{u—z)b{z)dz 

hJu h hj^ 

To deal with the first term we see that \h~^{a{u + h — z) — a{u — z))b{z)\ < 

Il'^'ll[-i5,<5] ll^ll[«-<5.«+i5]: so that the the first term converges to as /i — 0. As 

clearly the second term converges to a{0)b{u), we conclude that 
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lim - / a{u + h- z)b{z)dz = a(0)6(u) (16) 



Therefore, (jT4| - (jl6p show our result. □ 

Remark 2.4 The integrability condition for b in Lemma [K7\ is authomatically 
satisfied if b ^ C([0,oo)) 

Next result gives conditions under which m is differentiable, along with useful 
expressions for its derivatives. 

Proposition 2.1 Let m be the locally bounded solution of fO]) with f the den- 
sity function of an absolutely continuous nonnegative random variable Z with 
distribution function F. Assume that F and v are in C^([0,oo)) We have that 
m G Cg ([0, oo)) and for all u > 

m'{u) - ~ y)f'{y)dy + m{u)f{0)^ + v'{u), (17) 

m"{u) = (f>( [ m'{u - y)f'{y)dy + to(0)/'(m) + m'(u)/(0) ) + v"{u). (18) 







Moreover, 

to(0) = 'y(O), m'(0) = (/.u(0)/(0) + v'(0), and (19) 
m"(0) = </>(i;(0)/'(0)+m'(0)/(0))+t;"(0) (20) 
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In addition, if we define the functions 



wi{u) 



(j)m{0)f{u) +v'{u) 



(21) 



W2{u) 



(j)m'{{))f{u) + w'i{u). 



(22) 



We have 




■ u 



m! {u) — (f> 



m'iu - y)f{v)dy + wi{u), 



(23) 



■ u 



m"{u) = (p 



m'^u - y)f[y)dy + W2{u). 



(24) 



Proof. First of all, from ^ , the fact that m is locally bounded and the continuity 
of V we deduce that to is continuous on [0, oo). Also, the first equality in (|19p 
follows by ([5] Secondly, making the change of variable z — u — y in the integral 
contained in we can write 



Differentiating the previous expression (recall Lemma |2. II and Remark 12. 4p . we 
obtain 



this shows that m is continuoslly differentiable on (0, oo). If we make now the 
change of variable y — u — z \n the previous expression we obtain (jl7p . Also, 
taking limits in this expression as u J, 0, we obtain 



thus showing that to G Cq([0,oo)) and the second equality in ([T^ . Now, (|18p 
follows by differentiation in ([T7|) . and ([^0]) follows taking limits in (|18p as m J, 
Finally once it is shown that to G Cq([0,cx))), Lemma [2.11 and Remark 12.41 
allow us to make a straightforward differentiation in © to obtain ([25)1 and (IM)) . 
□ 





ni{z)f'{u — z)dz + ■m{u)f[Q) j + v'{u), u > 0, 



limTO'(M) = 0TO(O)/(O) + w'(0) = (/)u(0)/(0) + w'(0). 
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In order to obtain bounds of the derivatives of to, the following technical 
result will be useful. 

Proposition 2.2 Let TO2 be a function satisfying 

1712 {u) = (f) mi{u-y)fi{y)dy + vi{u), w > 0, (25) 

JQ 

in which < (f> < 1, mi and vi are continuous functions on [0, oo) and fi is a 
continuous function on (0, oo). Let 

/•oo 

Wi) ■■= / y'\fi\{y)dy, * = 0,1,2,... (26) 

Jo 

(a) ///o(/i) < 00, then 

||™2||[o,a=] < <?!'/o(/i)||mi||[o,:,] + 2;>0, (27) 

(b) If liifi) < 00, i = 0,1 we have, for all x > 

l|uTO2(u)l|[o,a:] < (-^0 (/l ) II ""T-i (u) 1 1 [o,:r] + ^1 (/l) 11^1 (^i) II [0,x]) + \\uVi{u)\\. (28) 

(c) If liifi) < oo, i = 0, 1, 2 then, for all x > 

||w^m2(M)||[o,j;] < {lo{fi)\\u'^mi{u)\\[o,x] + 2/i(/i)||utoi(m)||[o,:c] 

+/2(/i)||™i||[o,.]) + ll^iV(«)l|. (29) 



Proof. Part (a) is straightforward taking norms in (|25|) . To prove (b) and 
(c) we use that m" — {u — y + yY^ — J27=o ~ !^)"~*- Therefore, using 

([25]) . we can write 



(u-yY {u - y)y'fi {y)dy + (u) 



i=0 

from which we deduce easily 



|u"m(7/)||[o,,] < -^E ( J ll"""'"^i(")ll[o,.] 1 yViMdy + 

i=0 



Then, (b) and (c) follow easily from the previous expression applied to n = 1 
and n = 2, respectively. □ 

Taking into account the previous result and using the expressions (|23p and 
(j24p . next result gives bounds for the weighted derivatives of m. Later on (in 
Proposition l2.4p we will give sufhcient conditions in order to ensure the finiteness 
of these bounds. 

Proposition 2.3 Let m be the locally bounded solution to ^ with f the density 
function of an absolutely continuous nonnegative random variable Z , with distri- 
bution function F and finite variance. Assume that F and v are in Cq([0, oo)). 
We have the following 

(a) Let wi be as defined in (|2ip . We have 

Wm'W < 1^ (30) 
||W(.)|| < <^EZ\\m'\\ + \\uw,iu)\\ 

1-0 

11^2^/(^)11 < 4'i2EZ\\um'{u)\\+EZ^m'\\) + \\u^wi{u)\\ ^g^) 

l-cj) 

(b ) Let W2 be as defined in (|22p . We have 

ll-"ll < 1^ (33) 

||W(^.)|| < ^mm"\\ + \\uw,iu)\\ 

1-0 

11^2^.(^)11 < H2EZ\\um"{u)\\ + EZ^m"\\ + \\u^W2{u)\\ ^^^^ 



Proof. To show the bounds in (a), we use (|23p and apply Proposition [521 with 
m2 — mi — to', vi — wi and /i = / being a density function. Thus, /i(/) = 1 
and using (j27l) . we can write 

ll»7t'||[0,x] < 01 1 1 1 [0.x] + <X < OO, 
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As by Proposition l2 . 1 1 m' is a continuous function on [0, oo), then ||m'||[o.2;] < oo, 
from which we deduce 



(1 - (/))||m'||[o.a;] < ll^ill, Q<x<oo, 



Taking limits as a; — > oo in the previous expression, (1301) holds true, pip and 
are shown in a similar way, using and (I^Hl) . respectively, and taking 
into account that Ii{f ) — EZ''', i = 1,2. The proof of part (b) is similar, taking 
into account (I24p and applying Proposition 12. 2[ with TO2 = "ii = tiT'" , vi — W2 
and/i = / □ 

Our next task is to give conditions for m in order to ensure that m £ T>i. 
First of all, we state a technical lemma in order to simplify our hypothesis. 

Lemma 2.2 We have 

(a) Let vi e C([0, cxd)). // ||it^t)i(M)|| < oo, then \\vi\\ < oo and \\uvi{u)\\ < oo. 

(b) Let h e C((0, oo)). Let /^(/i) be as defined in < oo, i = 0,2, 
then Li (/i) < oo. 

Proof. To show part (a), let vi G C([0,oo). The continuity implies ||wi||[o.i] < 
oo and the result is immediate using the the following bound 

||u'wi(u)|| < lki||[o,i] + ll"^wi(")ll(i,oo) < oo i = 0, 1 
Part (b) is due to Cauchy-Schwartz's inequality, as using that = 



/i(/i)=/ u\h\{u)du<[ u^\h\{u)du] / \h\{u)du] = (/2(/i)/o(/i))'/' 



u\h\^/^{u)-\h\^/\u) 







□ 



Now we enunciate the main result of this section. 
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Proposition 2.4 Let m be the locally bounded solution to Q with f the den- 
sity function of an absolutely continuous nonnegative random variable Z with 
distribution function F. Assume that 

1. Z has finite variance. 

2. F and v are in Cq([0, cxd)). 

3. f andv" are in {[Q ., oo)) . 

4- Ii{f ) < oo, i = 0, 2, where ) is defined in (|26|) . 



5. \\u'^Wi(u)\\ < oo, i = 1,2, where wi, W2, are defined in (j2ip and (1221) . 
respectively. 

6. \\u'^w'l{u)\\ < oo, i = 1,2. 
Then we have: 

If Condition 2 is satisfied and in addition \\w2\\ < oo, then, 

||m"||<|^<^ (36) 
// Conditions 1-6 are satisfied, we have 

< 0((/o(/'O + l/'(o)l)h'™'H|| + 2/i(/'0l|W(u)i| + /2(r^ 

+ (t)f{Q)\\u^m"{u)\\ + \\u^w'l{u)\\ < oo (37) 
||^i'm-(u)|| < cb ((/o(/") + \f'm)\\u'm"{u)\\ + 2h{f")\\um" {u)\\ + /2(/")ll™"|| 
+ (f>f{0)\\u^m"'{u)\\ + \\u^w'^{u)\\ < oo (38) 



Proof. Note firstly that ([55]) is obvious by Condition 2 and Proposition 12.31 
(b) . Secondly, note that Condition 2 allow us to apply Proposition 12.11 Our 
starting point to prove ([57)) - ([55)) will be the expression for m" given in (ITO)) . 



12 



Moreover, as by condition 3 /' is differentiable, the integral appearing in this 
expression can be rewritten, by means of an integration by parts as 

m'{u~y)f'{y)dy^ ( m{u ~ y)f"{y)dy - m{0)f'{u) + miu)f'iO) (39) 
Jo 



Inserting (p9)) in (IT8| , we obtain 

m"{u)^<p( f m{u-y)f"{y)dy + m{u)f{0) + m\u)f{0)]+v"{u). (40) 







We differentiate the previous expression (note that by Proposition 12.11 m G 
Co([0, oo) and, by Conditions 3 and A, b — f" satisfies conditions of Lemma 
12. ip . Thus, we obtain for all u > 



2 



m"'(u) ^^yj m'{u~ y)f"{y)dy + m(0)/"(w) + ^ m^"^ {u)f^^-'^ (0) j + v'"{u) 

= <j) (^J\\u - y)f'\y)dy + J2m^'\^)f^'~'H0)^ +<("), (41) 
where wi is defined in (|21l) . Thus (|4ip verifies (PH)) in Proposition [521 with 

2 

m2 = m'", mi = m', /i = /" and ui = 0^mW(u)/(2-»)(O) + u;'/(u), 
SO that we deduce from Proposition [52] (c) 

\\u^fn"'{u)\\ < 0(/o(/")|k'm'(«)|| +2/i(/")||W(u)|| +/2(/")||m'||) 

2 

1=1 

and the first inequality in (|57)) follows easily from the previous bound. To show 
we use Lemma to differentiate (HI]), thus obtaining 



m™(u) = (^y^ m"(w - y)riy)dy + ™'(0)/"(w) + ^ m^'+D (w)/^^-^) (0) j + <(u) 
= [ /"m"(« - y)/"(y)ci2/ + ^ m^^+i^ (t.)/^^-') (0) ) + u;^'(u), (42) 



V i=l / 

where W2 is defined in (|22p . Then, as (fi5]l verifies pOI) in Proposition [52] with 

2 

7712 = mi = m", /i = /" and vi = </>E w^*^^Hw)/'^~'-'(0) + < ("), 

13 



we obtain from Proposition [52] (c) 

h'm™(u)|| < (/.(/o(/")|k'm"(w)|| +2/i(/")||W(u)|| +/2(/")||m"||) 

2 

thus showing the first inequahty in ((38)) To show the finiteness of p7l) and (p8)) . 
we win prove that 

Conditions 1,2 and 5 ^ ||w-'u'i(u) || < oo, i = 1,2, j = 0,1,2 

^ llu^rnW (m) II < oo, i = 1,2, j = 0,1,2. (43) 

To show let i = 1,2 be fixed. Condition 2 imphes that Wi G C[0, oo). 

Thus, by condition 5 and Lemma [2.21 (a) we have the first imphcation in (^5)) . 
For the second imphcation we use condition 1 and apply Proposition 12. 31 fa) for 
1 = 1, whereas for i = 2 we apply Proposition 12.31 (b). 
Now, note that 

Condition 4 ^ /,(/") < oo, i = 1, 2, 3 (44) 

which is immediate by Lemma [2.21 fb). Thus, using (|43l) . (HU and Condition 
6 we show the finiteness of ([57)) . Similarly, the finiteness of the bound in psp 
follows using ([^5|. ((3D), ((57)) and Condition 6. This completes the proof of 

Proposition [Ol □ 

As an immediate consequence of Proposition 12.41 and Theorem 12.11 we have 
the following. 

Corollary 2.1 Let m he the locally hounded solution to ([9]) with f the den- 
sity function of an absolutely continuous nonnegative random variahle Z with 
distribution function F. If conditions 1-6 in Provosition 12.41 are satisfied then 
m £ 2?!, withVi as defined in (1131) . Therefore, the approximation m, t > 0, 
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as defined in (jlip - (|12p verifies 




m-m\\< —\\m"{u)\\ + —\\um'"{u)\\ + 
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||M^m*''(u)|| < oo. 



16^2 



3 Application: Approximations for ruin proba- 
bilities 

In this section we will apply the results given in the previous one to ruin proba- 
bilities in the classical risk model, which are a well-known example for functions 
satisfying a defective renewal equation. First of all we recall how the classical 
risk model is defined (see [U Ch. 4] or [71 Ch.5.3], for instance). We con- 
sider an insurance company in which insurance claims follow a Poisson process 
{N{t), t > 0) with intensity A > 0. On the other hand, the individual claim 
amounts (Xi)igN* are identically distributed and positive random variables with 
finite mean, independent on {N{t), t > 0) . Suppose that the initial capital of 
the insurance company is U(0) := u > and it receives premiums at a constant 
rate c. In this setting, the probability of eventual ruin ^j{u) is the probability 
that the wealth of the company is ever negative, i.e. 



Call fi := EXi, and assume that /i > 0. The condition for no sure ruin is 



We will assume this condition from now on. Usually the ruin function cannot be 
evaluated in an explicit way. However, it is well-known that the ruin function 
satisfies a defective renewal equation (cf. 4, p. 105], for instance). In fact, let 
Fx be the distribution function of the claim amounts, and Fx := 1 — Fx- We 




(45) 




(46) 
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have 

^{u) = r ^{u - y)^^dy + r ^^dy] , > 0. (47) 

and therefore ^ holds true for the ruin probabihty, with 

F'iu) = f{u) = and viu) := ^ / f{y)dy = ^{1 - F{u)) (48) 

Note that / is a weU-defined density corresponding to the so-called equilibrium 
distribution of X (cf. [11, p. 14]). The rest of the section is divided in two 
parts. In the first one we will check that if the claim amounts are mixtures 
of gamma random variables with shape parameter a > 1 (and arbitrary scale 
parameter), then ^E* G so that Corollarv 12.11 holds true. In the second part, 

[21 . . 

we will give a method to compute Ml 5*, paying special attention to mixtures 
of gamma claim amounts. 

3.1 Conditions for optimal order of convergence in the ap- 
proximated ruin probability. Applications to mixtures 
of gamma claim amounts 

Our first result gives sufficient conditions in order that G 2?i . As an immediate 
consequence we will obtain sufficient conditions for mixtures of gamma claim 
amounts. 

Proposition 3.5 Consider the classical risk model, where the claim amounts 
X have distribution Fx- Assume that Fx S C'^[0, oo) and 

a) X has finite third moment 

b ) F^ verifies that lim Fx (u) exists and is finite. 

c) I,{F'^) = u'\F'j,{u)\du < ^, z = 0,2 

d) \\u^Fx{u)\\ < oo, and \\u^F^^\u)\\ < oo, i = 1,2,3 
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Then, the associated functions F and v for the ruin probability as given in 
(|48p satisfy conditions 1-6 in Provosition \2A\ 

Proof. We will consider the renewal equation for the ruin function as given 
in (|Tfl) and and check all the conditions in Proposition 12.41 
Condition 1. Let Z be the random variable whose distribution is F , as given 
in (|48p . To show that Z has finite variance, we note that an integration by 
parts shows us that EZ"^ = EX'^ / {2>EX) (cf [11, p. 15]) so, that if a) is satisfied, 
Condition 1 holds true. 

To check the rest of conditions, taking into account pS)) . the functions for 
the renewal equation and their respective derivatives are 

F [u) = f(u) = , and F [u) = f (it) = 

v{u) = ^il-F{u)), v'iu)^(f>^^^^ and v"{u) = cb^^^ (49) 

^J. n 

Condition 2. Taking into account ()49p it is immediate that if b) holds true, then 
F and v are in Co([0,cx))). 

Condition 3. Taking into account (j49|) . both F" and v" are in C^[0,oo), as 
F'x e C2[0,oo). 

Condition 4- We need to chack that hif") < oo, i = 0,2. This follows from c), 
as /" = ~F^/^i. 

To show Conditions 5-6 we recall (|2T]) and use (l49l) and (IT9l) to write 



wi{u) = 0w(O)/(u) + v'{u) = -0(1 - (j))^^^ (50) 

A* 

Now we recall (HU and use (gH) and to write 

W2iu) = H'{Q)f{u) + K(«) = 0V/(O)^^^ + <^(i - 4^)^^ (51) 

Recalling ([TOl) and using we have that 

(0) = 0/(O)w(O) + w (0) = 
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and using the previous equality, (j5ip and (1501) we can write 

= + 0(1-0)^. (52) 

Condition 5. We need to show that [[^^^^(u)!! < oo,i = 1,2. Taking into 
account d), we have that < oo, which imphes, recahing (|50|) . that 

||u^fi;i(w)|| < oo. Using the previous bound, ([52)1 and the fact that, by d), 
Iju^i^jf (m)|| < oo, we have that ||u^W2(u)|| < oo 

Condition 6. We need to show that [[^^^"(m)!! < oo, i = 1,2. This follows 
easily by ^ and d), as 

\\u'w'l{u)\\^^^^\\u^F^{u)\\<^ 

Ww'i{u)\\ < ^:^||.V/(u)|| + ^^}—^\\u'F^'iu)\\ < oo 
This completes the proof of Proposition l2.4l □ 

Corollary 3.2 Assume that, in the classical risk model, the claim amounts are 
mixtures of gamma random variables T{ai, Pi), with mixing weights {pi)"^i, that 
is 

n 

J^XH =^P^i^a.,ftH, (53) 
1=1 

where Fa^^jj^ are distributions having density fai.pi defined in ^ , pi > and 
Pi + ---+Pn = l. Assume that ai > 1, i — l,...,n. Then, the ruin probability 
satisfies that Tp G Vi, with 2?i as defined in (|13p . Therefore, the approximation 
mI'^^iP, t > 0, as defined in (|lip - (jl2l) has uniform order of convergence 1/t^, 
that is 

IIMfV - V'll < ^\W'{n)\\ + ^||W/"(^^)II + Y^lkV"(^^)ll < ^. 



Proof. To prove the result we will show that Fx satisfies the conditions 
of Proposition 13.51 Condition a) is obviously satified. To check the rest of 
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conditions, we will use the following simplification. Denote by Fa := F^^i, a 
gamma distribution T{a, (3 = 1). Recall that a gamma distribution r(a, /3), has 
distribution function Fa^piu) = Fa{l3u). Therefore, we can write (j53p as 

n n 

Fx{u) = ^piFa,{l3iu), and Fx{u) = ^p,Fa^{l3,u) 

and it is clear that if F^, a > 1 satisfies conditions b)-d) in Proposition l3.5[ then 
Fx will also do. Then, we will check conditions b)-d) for Fq, with a>l. Note 
firstly that Fa is infinitely differentiable on (0,cxd). In particular, its density 

= f^e'"""~\ «>0. (54) 

satisfies condition b) in Proposition [XU whenever a > 1. For the rest of condi- 
tions, note that 

F'J,{u) = -^e-^u'^-\a-l-u), w > 0, (55) 

F'^'{u) = -^e-^u''-'^{{a-l){a-2)-2{a-l)u + u''), w > 0. (56) 
r(a) 



Observe that the previous equalities follow by differentiation in ([54| for a > 1, 
and are still valid for a 1 . To show condition c) in Proposition 13. 5[ note that 



r(a) Jo r(a) Jo 







T{a) T{a) 
To show Condition d), note that it is clear from ([54)) -([56 | that 

lim u'^F^'^ (u) = 0, i 1, 2, 3 (57) 

and using F'Hopital'rule applied to v?'Fa{u), we can write 

lim u^Faiu) = lim 2u^F^{u) = (58) 



From ([57|) . (f58| and the fact that Fa and are continuous and bounded at 
the origin we obtain that |jit^FQ(u)|| < oo and ||u^F^(m)|| < oo. To show that 
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< oo, i — 2,5, we take into account ((57)) and note that 

lim v^F'^{u) =0, a > 1, 

lim ^^^-N, ^ ~ lim = 0, a>l. 

Note that the last equality follows as, for a = 1, a — 1 = 0, whereas for a > 1, 
lim„_i.o = 0. Thus conditions a)-d) in Proposition 12 .41 are verified, and the 
conclusion follows by Corollary 12.11 □ 

3.2 Numerical computation of the approximated ruin prob- 
ability. Applications to mixtures of gamma claim amounts 

In this Section we will give a method to compute the approximated ruin proba- 
bility, using that non-ruin probability is the distribution function of a geometric 
sum. Our approach is based on the following representation of L^g, as defined 
in ([T]), when g :~ Fx is the distribution function of a nonnegative random vari- 
able X. In this case, the approximation L^Fx can be rewritten in the following 
terms (cf. P]). Let (t>x{- ) be the Laplace-Stieltjes transform of X , that is 

^x{t) := Ee''^ = [ e-*"dFx(u), t > 0. 

J[0,oo) 

We define a random variable X'* taking values on k/t, fc € N, and such that 
P(X'' = k/t)^tI}^^P{t), fceN, (59) 



Let {S{t), i > 0) be a collection of gamma random variables r{t, 1), as consid- 
ered in the Introduction. The following equality holds true (see [I]), 

P{X-<u)JY.^-^<t>'x\t)-EFx(^^^^±^), u>0, (60) 

k=0 ' ^ ^ 

As mentioned in the Introduction, the first equality is the so-called Widder's 
formula in [5]. Let L^Fx be the approximation defined in Recalling Q we 
see therefore that 

/ j-\[t«] + l ^lUu]) 

^*^^^"^^r(RTT)^^ W = W*<^^) (61) 
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Thus, L^F can be obtained either by straightforward differentiation using the 
first equality or by computing the probabihty mass function of X** (second 
equality). Now, we will see the computational advantages of using the second 
method to approximate -0, the ruin probability in the classical risk model, as 
defined in Section 3.1. To this end, denote by the non-ruin probability, that 
is '0 = 1 — ■(/'. It is well-known (cf. [4l p. 104], for instance) that 

M 

4>{u) = L, < u) (62) 

1=1 

in which M is a geometric random variable, with probability of 'succes' p = l—(f>, 
where is as in 1^. That is, P{M ^ n) = _p(l-p)" = (1-0)0", n = 0, 1, . . . , 
and (Li)igN-> is a sequence of i.i.d. random variables, having density / as defined 



in (|48p . and independent of M. Applying (|6ip to non-ruin probability and 

/ M \ •* 

taking into account that [ Li ] has the same distribution as J^fLi (see 
[SI Proposition 2.1.]) we can write 

Lt^iu)=P\(Y,Li] <u] =p(j2Lf , ">0 (63) 



\i=l / I \i=l 



and using the accelerated approximation M\ ip as defined in (jll[) we can write 

Therefore, we can approximate the ruin probability by evaluating the distribu- 
tion function of a discrete compound geometric distribution. This allows us to 
use well-known evaluating techniques for compound discrete distributions (Pan- 
jer's recursion, for instance, see [4| p. 50])). Approximations for ruin probabilities 
by means of the discretization of the summands in (j62l) have been proposed in 
the literature (cf. 4, p. 110]). Perhaps the most natural way to discretize a 
random variable is to round it from below or from above. However, rounding 
methods are difficult to apply when the distribution function of a random vari- 
able cannot be given in an explicit way (consider a gamma random variable 
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with a shape parameter being not a natural number, or its equihbrium distribu- 
tion, for instance). The computational advantage of our method is that we can 
evaluate the probability mass function of Lf, whenever the Laplace-Stieltjes 
transform of the claim amounts is known. The expression for the discretized 
record lows and their behaviour when dealing with mixtures, are collected in 
the following. 

Proposition 3.6 Consider a non-negative random variable X with distribution 
function Fx and Laplace-Stieltjes transform ^x- Assume that X has finite 
mean fi. Let L he a random variable having the equilibrium distribution of X , 
that is whose density is given as 

fdu)^^^, u>0 (64) 
M 

We have the following. 

(a) Let L** the discretization given in ([5^. We have 




(b) Assume that Fx is a mixure of random variables 

Fx =PiFi + ---+PnFn 

where the mixing distribution functions (^i)"=i; have finite mean 
and {pi)^^i are the mixing probabilities. 

Let (ii)"^-^ be random variables having the equilibrium distribution of{Fi)"^^. 
Then L** as given in (I59p verifies 

p^^.t ^k^^ PH^p^Ll' = ^) + . . . + 2ll!^P{L',^ ^ ^) (65) 
< /i t /i t 
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Proof. To show (a), we use an integration by parts to write 
<fL(t) = / e-'"^^du - (1 - <5x(0) 

J[0,oo) M 

and therefore, applying Leibnitz's differentiation rule, we can write 
Thus, using ()59p . we can write 



thus proving (a). To show (b), note that we can write 
Fx_ _ 'PUhF\ ^ PiUi K 

and, taking into account the previous expression we have 

Jlo.oo) MM M 



Thus, 



= Pl!^P(^L[' = -) + ■■■ + ?^p{m 

j.1 t /it 

thus showing part (b) □ 

As an immediate application of the previous result, to mixtures of gamma 
random variables, we have the following. 

Corollary 3.3 Consider X a nonnegative random variable, let L be its equilib- 
rium distribution as given in (|64[) and consider its discretization L** as given in 
(|59p . We have the following 



(a) Assume that X is a gamma random variable T{a,l3), that is, having density 
as given in ^ . Consider the cumulative distribution of a negative binomial 
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random variable with a > 'successes' and probability of 'success' p, that is 
CDF.NB{k;a,p) =^ ('^ ^ "^Vl - fc = 0,l,--- 

Then, 

PiL'^^'-)^l(l-CDF.NB(k;a,p=^)) (66) 

(b) Assume that X is a mixture of n gamma random variables, with mixing 
weights i^at is 

Fx =PlFi + ---+PnFn 

in which each Fi has distribution T(ai,Pi). Then, 

f^p, A - CDF.NB (k; a,, 





P(L-* = -) = — 7 r — (67) 

t / piai pna„ 



t 



Pi Pn 



Proof. To show (a) note that the Laplace-Stieltjes transform of a T{a, (5) 
random variable is 



and therefore, 



t + P 



Recalhng that p — a/P and applying Proposition 13.61 (a) . we have 

which shows (a). Part (b) is immediate by part (a) and Proposition 13.61 (b). 
taking into accomit that pi — ai/Pi. □ 

Remark 3.5 As shown in [T], when X is a gamma random variable T{a,P), 
the weights of X'* , as defined in (|59p correspond to the ones of a negative 
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binomial random variable. From Corollarv \'6.'S\ (a) we deduce that the discretized 
equilibrium distribution of a gamma T{a^ (3) random variable is constructed by 
cumulative sum of the afore-mentioned weights. 

Example 3.1. Approximation of ruin probabilities when the claim 
amounts are mixtures of gamma distributions In this example we show 
some numerical approximations of the ruin probabilites with the method de- 
scribed above, by considering mixtures of gamma claim amounts. First of all 
we describe the steps needed to build the approximation 

1. Computation of P{L'* = j) by the mixture formula given in (|67l) . For fixed 
t and fixed values of (pi)-Lij (Q!i)"=i and (/3j)"=u '^^ need the probability 
distribution of the corresponding negative binomials. In our case, we used 
MATLAB to generate these values, for k — 0,1,2, . . . 



2. Computation of L^ip, using (|63)) . Note that this can be done using Panjer's 
recursion, a popular method for evaluating compound distributions. Panjer's 
recursion applied to the geometric sum given in (|63p provides the following 

M 

recursive formula for evaluating the probability mass function of L** (cf. 
12 P- 50]) 

1-6 



i=l 



1 - <j)P{Lf = 0) ' 



if fc = 0: 



M 

t 




1 - (j)P{L'* = 0) 

In our case we used an EXCEL worksheet to generate, for fixed t, and 
the computations given in the previous step, the corresponding values of 
L*Mklt). 



if fc= 1,2, 



25 



3. The final approximation for the non-ruin probabihty is 

Note that this implies to repeat steps 1 and 2 for t (second term above) and 
2t (first term above) 

If the shape parameters {ai)f^i of the initial mixture of claim amounts are 
chosen to be a; > 1, i = 1, . . . ,n, CoroUarv 13.21 ensures us a uniform order of 
convergence of . 

We give three numerical computational examples. For them we use (f> = 0.9 
and t — 5. Results are shown in Table [T] 

• First column provides us the approximation of non-ruin probability for 
exponential claim amounts, having mean 1 that is, T(a — 1, jS = 1). This 
example can be used as a test, as exact non-ruin probabilities can be 
computed in an exact way by the formula (cf. 4, p. 93]) 

Vi(M) = 1 - (1 -p)e-P", u>0, 

where p — 1 — (f>, as above. Note that this function was the used in :9, Ex- 
ample 2.1] for numerical computations (recall Remark l2.3l) . and provides, 
with the given parameters, exact values up to four decimal places. With 
respect to Column 1 in Table 1 the approximation is recalculated using 
the Panjer's recursion method described above, whereas in [5] the explicit 

[21 ~ 

expression for Mg -0 ^^-s used. 

• Second column provides us non-ruin probability for Gamma distributed 

3 

claim amounts T{a ~ 2'^^' "^^^ interest of using the approximation 
in this case is that, when a ^ N, there is no explicit expression for the 
ruin probability. However, alternative approximate expressions can be 
obtained by series expansions (see \TU\ Y 
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• Third column provides us non-ruin probability for claim amounts being 
a mixture of the previous cases, with weights pi = P2 = 1/2. We have 
chosen the same scale parameter in both terms, to easen comparability, 
but note that there is no computational problem in choosing different scale 
parameters. 



u 



.J 
-I 

20= — 

150 



Exponential claims 



Gamma a = 3/2 claims 



Mixture 



30= 



40= 



2D0 
_5_ 



M^'^V(l) = 0.1856 
A^^V^S) = 0.4538 
M^^^ipilO) = 0.6677 
Aff lVi(15) = 0.7975 
M^^lV;(20) = 0.8766 
Ml^^tp{30) = 0.9553 
M^tpiAO) = 0.9854 



M, 



M, 



M. 



M, 



^{l) = 0.1648 
1^/5(5) = 0.3940 
''^^(lO) = 0.5949 
''^5(15) = 0.7248 
'^i){2Q) = 0.8190 
'^'4){2,Q) = 0.9191 
■lV;(40) = 0.9639 



Mf 1^5(1) = 0.1726 
Mf^{b) = 0.4159 
Mf^{lQ) = 0.6225 



m|^V(15) 
M]^lt^(20) 



0.7560 
0.8423 



V(30) = 0.9341 
^^(40) = 0.9725 



Table 1: Approximation of non-ruin probability for different claim amounts 
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